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Abstract 

We consider a linear finite spring mass system which is perturbed by 
modifying one mass and adding one spring. From knowledge of the 
natural frequencies of the original and the perturbed systems we study 
when masses and springs can be reconstructed. This is a problem about 
rank two or rank three type perturbations of finite Jacobi matrices 
where we are able to describe quite explicitly the associated Green's 
functions. We give necessary and sufficient conditions for two given 
sets of points to be eigenvalues of the original and modified system 
respectively. 
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1 Introduction 



We study a problem on inverse spectral analysis of Jacobi matrices. Our mo- 
tivation is to understand the behavior of oscillating Micro Electro Mechan- 
ical Systems (MEMS) which are often modeled as masses mo, mi, ...,mN-i 
joined by springs with stiffness parameters (elastic coefficients) ko, ki, /cat 
and equilibrium lengths lo,li, ...Jn- 
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Figure 1: mass-spring system with 3 masses 

The first and last springs could be attached to fixed walls as in Figure [H 
or free, in which case the model is the same just setting the elastic coefficient 
of the spring at the free end equal to zero. The masses are allowed to move in 
the X— horizontal direction, with no friction and in absence of external forces. 
Using Hooke's and Newton's second laws, after normalization the following 
equation is obtained: 

where the entry z/j(t) of the iV- dimensional vector z7(t) describes the position 
of mass rrii at time t (see [1],|5]) and J is the Jacobi matrix 
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from now on we 
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parameters 7j = of the springs instead of their length li and the Hooke's 
coefficients k^). In the last years several experimental papers [n],Iin] were 
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written on the possible methods of determining micromasses with the help 
of oscillating microcantilevers, using the spring-mass system approach. A 
possible theoretical basis for such a task is given by the paper of Y.M.Ram 
^ (1993), who considered the inverse spectral problem of reconstructing the 
Jacobi matrix ([1]) by its spectrum and the spectrum of the perturbed matrix 
J with 



Oat-I — CLN-1 + — • 

tun-i mN-i 



JN-2 



rfiN-i 



JN-2-, 



rriN-i > 0, 7 G 



(3) 

the other entries remaining without change. He obtained the necessary and 
sufficient conditions for two point sets to be the spectra of such a pair of 
matrices J, J and provided a method of reconstructing the matrices by the 
spectral data. Ram's results were partially extended by P. Nylen and P. 
Uhlig in [n],[T] who considered the case of an analogous interior perturbation 
affecting the entries a„, Namely, they study the problem of changing 

the mass m„ by m„ for a fixed n, 0<n<A^ — 1, and attaching to it a new 
spring of elasticity parameter k, the other end of the spring being fixed: 



ko 



A k 
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Figure 2: Perturbed mass-spring system 
In this case the coefficients a„, 6„_i, &„ are modified as follows: 

fin = -r-(7n + 7n+l + l) = O"^ {cLn + M), In-l = 0bn-l, K = dhn, (4) 

where 

e = e{n):=.f^, M:=^, 7 := y (5) 

V m„ Win I 

rhn denotes the perturbed mass and I the length of new spring. All other 
entries of J remain unaltered. 

The inverse problem for such a perturbation is the problem of recon- 
structing J when we have the eigenvalues Ai < A2 < . . . < A^r and Ai < A2 < 
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. . . < Aat of J and J, respectively, where J = J{n) is the perturbed matrix 
with the changed coefficients (jl]) and all other entries the same as 

P.Nylen and F.Uhlig obtained in [6] necessary conditions for the spectra 
of the matrices J and J and offered a method for reconstructing the possible 
matrices in the cases where there is a finite number of solutions. However, 
as it will be seen later on, their conditions are not sufficient and the recon- 
struction does not give all the possible matrices. 

The goal of this paper is to develop the direct and inverse spectral the- 
ory for the interior perturbations of the described type (jl]), to obtain the 
necessary and sufficient conditions for the spectral data and the complete 
description of the possible Jacobi matrices with such spectral data, provid- 
ing an explicit algorithm of reconstruction. 

The main results are theorems [3] (necessary and sufficient conditions) and 
m (description of the possible matrices). The algorithm of reconstruction is 
given in their proof. In forthcoming papers we plan to give the complete 
solution of this inverse problem for seminfinite matrices. 

We shall use the following notations 

Am„ = m„ - m„, K := — ^ (6) 

Arrin 

and assume that Am„ > 0, i.e. the perturbed mass m„ is greater than the 
initial mass m„. 

We shall have the spectra of the initial and the perturbed Jacobi matri- 
ces cr(J) = {Ai, A2, Atv}, cr(</) = {Ai, A2, Atv}, a given integer number 
n G {0, . . . , N — 1} which indicates the place the mass is modified and the 
parameter of the perturbation K = as the spectral data for the inverse 

problem of reconstructing J and J. Notice that multiplying all the masses rrii 
and elasticity parameters 7j in ([2]) by the same constant does not change the 
Jacobi matrix ([1]) and the frequency characteristics of the free oscillations, so 
we obtain a mechanically "equivalent" system. Thus, we cannot reconstruct 
exactly the masses and the elasticity parameters from the frequencies. How- 
ever, their fractions ^, ^ = 9'^ , ^ = M will be found. That is the 
reason why we will only work with fractional parameters of perturbation 6"^, 
M, K instead of the masses mi,mn and elasticity parameters 7^,7. 

^When n — we have only the two last equalities in (jlj since &„_i is not defined. 
Analogously if n = — 1 we have only the first and last equalities in @ , since 6„ is not 
defined. 
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Taking into account and we immediately have 

M = {0-'^ - 1)K, (7) 

thus we always know one of the three parameters of perturbation 6'^, M, K 
from the two others. 

We should remark that there is a huge variety of inverse spectral problems 
and many of them have been deeply studied. We refer the reader to the books 
[2] and |1] for valuable information on this important subject. 

2 Preliminaries 

The eigenvector c of the matrix J, which satisfies 

Jc= Ac, (8) 

can be written after normalization as c = (Po(A), -Pi(A), -P2(A), . . . , P/v_i(A))^ 
with Po(A) = 1 where -Pi(A) is a polynomial in A of degree z G {0, 1, . . . , — 
1}. From ([8]) we get the expressions: 

ao + 6oi^i(A) = A, 
6i_ii^_i(A) + aiPi{\) + hPi+i = APi(A), (9) 

hN^2PN-2{\) + aN-lPN-l{\) = APiV-l(A). 

Therefore the polynomials Pj are defined by the conditions 

Pi+i = Ui^-ai)Pi-bi_,Pi_,}, P_i = 0, Po = l, (10) 

i = 0, N — 2. It follows from the last equation in (j9]), that if A is an 
eigenvalue then Qn{^) = 0, where we define 

g^(A) := AP;v-i - aiv-iPiv-i(A) - 6jv-2P7v-2(A), (11) 

and each root of the A^— degree polynomial Qa? is an eigenvalue of J, too. 
The polynomial Qat is equal to the characteristic polynomial of J times a 
constant. 
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Let Pi and he defined as in ([H]) , ( ITT]) and let Pi, Qn be the corre- 
sponding polynomials for the perturbed operator J. We now try to get an 



expression that relates Qn and Qn- 

Fix n G {0, 1, . . . , — 1} and define the polynomials (pi for n < i < as 
follows: 

1 1 

V^n = (^„_^i = — , ^i = 7 {(A - ai_i)(y9i_i - 6j_2</'i-2} , (12) 

On k-l 

where we set h^-i = 1- Note that the degree of the polynomial yjj is i — n — 1. 
These are called polynomials of second kind while Pi are called of first kind, 
see [1]. 

Lemma 1. With the definition of (fi given above we have for any n G 



where A = A(^~2 - 1) - M and r{n) = e{n) (see ^) if n e {0, N - 1} and 
T{n) = 1 otherwise. 

Proof, a) Case < n < — 1 . First we shall prove 



forn + 1 < 2 < - 1. 

Using dl]) and the recurrence relations (Q for Pj and Pj we obtain 



{0,l,...,Ar-l} 




(13) 



Pi = ifiiAPn + Pi 



(14) 
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To obtain the last equation we added and subtracted AP„ — -^Pn- Therefore 

= ^AP^ + P„+i (15) 

= ifn+lAPn + Pn+1, (16) 

SO f lT^ holds with i = n + 1. 
For i = n + 2 we have 

Pn+2 = 7 |(^ ~ (^n+l)Pn+l ~ &ra-Pnl 



SO ( !T^ holds in this case. Assume (fT4|) holds for i — 2 and i — 1 where 
z — 2 > n + 1 and let us prove it holds for i. First, 



Pi = -^[{\- a,_i)Pi_i - h^2Pi-2] 



see (j9]), where the coefficients aj_i,6j_.i and 6^-2 are unperturbed since i > 
n + 3. Therefore 

P^ = -r^{{X- ai){^i_iAPn + Pi-l) - hi-2{^i-2APi + Pi_2)} 
= ^iAPn + Pi 

and (fT^ holds. To prove (fT3|) recall that 

QtV = (-^ ~ OAf_i)P/v_i — feAr-2-P/V-2 

= (A - aAf_i)((y9Ar_iy4P„ + Pn-i) " 6Ar-2 (v5iV-2^Pn + PiV-2) 

= (^nAPu + Qjv 

and (IT^ is proven. 

b) Case tt, = 0. Note that Pq = -Po = 1- Using the first equality in ([9]) 

Pi = i(A - ao) = {A - ^2(ao + M)} 

= ^^{P1 + ^1A}. 

Analogously we have 

P2 = ^ {P2 + V2A} 
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Using induction as in case a) we obtain Pi = 6 {Pi + fiA} for i = N — 1, N — 2 
and then Qn = ^(Qn + Aip nPqJ ■ 

c) Case n = N — 1. From (Q we know that 

-P/V-l = 7 {{^ " <^N~2)Pn-2 — bN-sPN-s} 

ON-2 

and 

Pn-1 = ~ {(A — aN-2)PN-2 — ^Af-s-P/v-s} = r Pn~i = '^Pn-i 

Therefore 

Pn-1 = ^Pn-i (17) 
Analogously and considering (fTTj) . we find that 

Q^(A) = (A - a^_i)P^_i(A) - bN-2PN-2{\) 

= (A - ^2(a^_i + M))^P^_i - ebN-2PN^2 

Adding and substracting 9\Pn~-i — ^Pn-i we get 

= OQNiX) + e{X{^ - 1) - m)Pm-i{X) 

Therefore 

Qn = e(^QN + A(PnPn-?) 
Observe that (p^ = 1 if n = N — 1. □ 
Let us define the jj Green's function 

where (.) denotes scalar product and 6j{i) = 1 if i = j and if i 7^ j, 

t,j e{o,i,...,N-i}. 

Lemma 2. For < n < N — 1 

LpN{z)Pn{z) 



G{z, n, n) 



Qn{.z) 



Proof. From the definition of polynomials P„ see or fllUI) we know that, 
for n > 1, Pn{z) = if and only if z is an eigenvalue of J[o.n_i] the n x n 
upper left corner of J, that is: 





( flo 




■■ 


. ...\ 




ho 


ai 


h ■■ 





J[0,n-1] — 
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■ ■ 


■ Ctn-l / 



Inductively from the definition it follows that 

1 



Pn(A) 



hobi ...h 



A" + lower degree in A 



n-l 



Since 6o^i • • • &n-i-Pn(A) and det(A — J[o,ri-i]) are monic polynomials of degree 
n and with the same zeros, they are equal, i.e. 



det(A - J[o,n-i]' 



bohi . . . bn^i 

Analogously from ( 1121) it follows inductively that 

1 



18) 



V3Ar(A) 



-A 



+ lower degree in A 



hN-2 ■ ■ - hn 

(recall &Ar_i := 1). If we define the matrix , n < N — 2, 

( On+l hn+1 \ 

hn+l hn+2 ' " ' 

hn+2 0,n+3 ' ' ' 



Jl 



ln+l,N-l] 



\ 







■■■ a^-ij 



then ipi is related to the matrix J[„+i jv_i] in the same way as Pj, defined 
by ([n]), is related to the matrix J. In fact the yjj are the Pi for the matrix 
J[n+i,N~i], multiplied by ^. 
We get similarly as (fTHjl 



V3Ar(A) 



1 



hn ■ ■ ■ hN-2 



det(A - J[„+i,Ar_i]) 



(19) 
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and 

Qjv(A) = 7 ^7 det(A - J) (20) 

Oq . . . 0^-2 



Therefore, if 1 < n < iV - 2, 

(/'Af(A)P„(A) _ det(A - J[o,n-i]) det(A - J[n+i,N-i\] 



Qn{\) det(A - J) 

When = 

V9jv(A)Po(A) _ det(A - J[i,Ar-i]) 

Qn{X) ~ det(A - J) 

and for n = N — 1 

V9Ar(A)PAr„i(A) _ det(A - J[0,Ar-2]) 

Qn{\) ~ det(A-J) ' 



(21) 



(22) 



(23) 



Recall (^iv(A) = 1 if n = - 1. 

According to Cramer's rule the solution of {z — J)u = 6n is the vector 
u = (m(0), u(1), ...,u{N - l)f with 

where {z — J)j is the matrix z — J with the j-column substituted by Since 
u{n) = {6n, {z - J)'^6n) = -G{z,n,n), we get 

Now observe det(2 — J)„ = det(A — J[o,„--i]) det(A — J[„+i_7v-i]) when 1 < 
n < N — 2. This can be checked expanding the determinant on the left 
side by cofactors with respect to the n-th column. When = 0, — 1 we 
get det(2r — J)o = det(A — J[i,Ar_i]) and det(2; — J)n~i = det(A — Jio,n-2]), 
respectively. Therefore from formulas fl?I]) to ([23]) the theorem follows. □ 

With the help of the two lemmas proved above, now we are ready to give 
an explicit formula for the Green function at a point nn. 

Theorem 1. 

rn ^ W 1 ^ nf=i(A-A,) ^^ 
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Proof. From lemmas [T] and [2] we get 

T{n)QN - Qn 



T{n)QN{X{e-^-l)-M) 
Similar to fl20l) we have 



G(A,n,n) (26) 



Qn{X) = ^ \ det(A - J) (27) 

Oo . . . 0Ar_2 

Using (jlj) and writing the determinant as a product involving the eigenvalues 
of J we get 



(r(n) defined in lemma [T]), and from (120 



1 ^ 

Qn = 7 ^ - (29) 



6o • • • &Ar-2 . , 
J = l 

Using (EH]) and ([29]) in (EE]) we get for < n < - 1 



rrin / nj=i(-^ ~ 



where we made use of ([5]). Multiplying both sides of the last equation by 
= 6'^ -1 and taking into account (^ we get (125]) . 

□ 

Remark 1. Multiplying both sides of fl30]) by A(Am„) — k we obtain, if 
is not a pole of (^(A, n, ra), 

"^n = ?7i„Q^f4^ "4 A = — ^ or G'(A,n,ra) = (31) 

n;=i(A-A,) Am„ 

From Lemma [2] and equations f[T5]) . flTIJl) we know that the roots of (^(A, n, n) 
are the eigenvalues of J[o,n-i] and J[„+i^7v-i]-It will be seen in Lemma E] that 
common eigenvalues of J and J are roots of G{X, n, n) too. So if we know one 
of these points, or the value of if not a pole, (for example not eigenvalue 

of J), plus m„ and the eigenvalues of J and J, then ihn can be determined 
from ( 131]) . A particular case of (|3T]) appears in |8], formula 2.17. 
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3 Direct problem 

Let 4'i,4'2, ■■■,i'N be an orthonormal system of eigenvectors for J with cor- 
responding eigenvalues Ai,A2,...,AAr. For example, one way to get a system 
of orthonormal eigenvectors of J is to consider the polynomials P„ defined 
by (fTOl) and normalize the eigenvectors {1, Pi{Xk), P2{^k), PN-i{^k)) cor- 
responding to eigenvalue A^. We obtain in this case lipkii)]"^ ~ — 

Lemma 3. 



N 



k=l 



2 



\Mj)\ 

\k — z 

Proof. Let 5j = Ylk=i(^ki^k, then ak = {Sj,ipk) = i^k{i) and 

^ ^ ^ 1 " \ib i 'M' 
(S^ akipk,{J - zy^y^^ai'ipi) = y]afcai(V'fc,T ipi) = 

k=l 1=1 k,l=l ' k=l 

□ 

The following Lemma can be proven using Lemma [31 The spectrum of 
the operator T will be denoted by cr(T). 

Lemma 4. Let /3 G o{J) = {Ai, A2, A„} 

a) IfO<n<N — 1, then G{f3, n,n) = 00 or G{f3, n,n) = 

b) If n = or n = N — 1, then G{(3, n,n) = 00. 
Proof, a) From Lemma [3l we know 



G{z,n,n) = ^ 



Xk — z 

k=l 

Assume (3 = Xk^ G cr(J). There are two possibilities: either ipkoin) 7^ 
or ipkoin) = 0. If the first holds, then Afc^ is a pole of G{z,n,n). In fact, 
taking the left and right limits along the real axis we get G^Xk^—, n, n) = +00 
and G{Xko+,n,n) = —00 respectively. If the second possibility holds, then 
the eigenvector ipko^ which corresponds to the eigenvalue of (3 = A^q of J, 
vanishes at n. This implies that A^q is eigenvalue of J[o^n-i] and Jin+i,N~i] 
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too. Therefore ipj^^X^g) = Pni^ko) = and using Theorem [2] follows that 
Gi(3 = \ko,n,n) = 0. 

b) Observe that Pq = I and similarly Pn-i 0, since otherwise QNi^ko) = 
P]S[-i{Xko) = and this would imply P, = for all i. Therefore, using 
\Mr)\^= E,Viy(C)P '^^°^''^ ^^'^"^ = ^ or N -Is^nd then G{/3,n,n) = oo 



12 _ |Pr(Afc)P 

follows. □ 



Using Lemma [3] and formula fl2^ we obtain the next formula, which will 
be used in the following results: 



1=1 

where 



N{X) := i^^^fi^ (33) 

n7=i(A-A,) 

and ipi,ilj2, ■■■,'ipN is an orthonormal system of eigenvectors of J with corre- 
sponding eigenvalues Ai, A2, \n- 

We call an eigenvalue A^ G cr^J) unmovable if it does not change after the 
perturbation, i.e. if A^ G cr{J) fl o-{J). 

Lemma 5. (Properties for the unmovable eigenvalues and 6) 

a{J) n a{J) = a{J) n (^{A : ^(A, n, n) = 0} U {K}^ (34) 

#|ct( J) n a{J) n {A : G(A, n, n) = 0}| < min(n, - n - 1). (35) 
A, G (t( J) n {A : G(A, n, n) = 0} ^ Ar(A,) = 6^ (36) 

A, = G (t(j) ^ N{Xr) = e^ + {i- e^)\Mn)\^ > (37) 

K ^ a{J) =^ N{K) = (38) 
K e a{J)Ua{J) ^ K e a{J)na{J). (39) 
If 1^ K , then \n o,nd if \i 7^ K, then Ai 7^ Ai. 
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Proof. "D" for (El. If A G a{J) but A ^ a{J) then A is a pole of N, by the 
definition of A^. 

If G(A, n, n) = 0, then from (|25|) . (|32|) . we get 



and A is not a pole of A^. Therefore, if one eigenvalue Xj^ of J coincides with 
one root of G{X,n,n), then some eigenvalue of J has to coincide with Xj^. 
Since the spectra of J and J are simple, no more eigenvalues coincide at that 
point. In case the eigenvalue Xj^ of J coincides with K, then from ([32D we 
get 



and Ajo = K is not a pole of N. Then some eigenvalue of J has to coincide 
with AjQ. 

"C" for If A G a{J) n (t(J) then A is not a pole of A^, by the 

definition of A^. It is enough to consider the case X ^ K. Using Lemma H] 
we know that A is either a pole or a zero of G{-,n,n). From fl5^ . poles of G 
that are not K, are poles of A^. Therefore A is a root of G if A G cr( J) fi o"( J) 
and X ^ K. 

Further, as mentioned in Lemma H] a), if an eigenvalue of J is zero of 
G{X,n,n), then it is a common eigenvalue of J[o,n-i] and J[„+i,Ar_i]. Since, 
according to |3j, there are at most min(?T,, A^ — 1 — n) of them, we conclude 
that the common eigenvalues of J and J are at most min(n, N — 1 — n) plus 
possibly the point i^', so f l35|) is proven. Notice, that by b) of Lemma HI 
cr(J) n {A : G(A,n,n) = 0} = when n = or A^ — 1, so the only possible 
common eigenvalue for J and J in this case is K. 

Implications f l5B]) -f l5I?]) are obtained using the proved part of the Lemma 
and formula fl5^ . 

Let us now pass to the last assertion. If the matrix A[o,Ar_2] is constructed 
from a Jacobi matrix A[o,Ar_i] by deleting the last column and row ( similarly 
first column and row), then cT(y4[o,Ar_2]) C {Pi, [3]^) where f3i, (3^ are, respec- 
tively, the smallest and largest eigenvalues of A[o,Ar_i]. See Corollary 2.5 in 
[3], for example. Using this fact with matrices J = J[o,Ar-i], J[o,n-2], ■■■ and so 
on, we get that the eigenvalues Ai and Aat of J, cannot be eigenvalues of the 
submatrices J[o,n-i], J[n+i,N-i] defined in Theorem [2|, and therefore cannot 
be roots of G{-,n,n). From Lemma [5] follows Aj. 7^ A,., r = 1, A^, unless equal 
to K. □ 



9^ = N{X) 



(40) 




(41) 
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Theorem 2. (Interlacing properties for two spectra) Let cr{J) = X2, ^n}, 
(t(J) = {Ai, A2, Aat}. Take p G {0, 1,...,A^} such that \p < K < \p+i, 
where we define Xq = —00 and Xn+i = 00, being K the parameter of the 
perturbation (jl]), ([6]). Then there is exactly one eigenvalue of J in each of 
the following intervals: 

[Aj,Aj+i), j = l,...,p-l, 
(Aj,Aj+i], j =p + 1,...,A^- 1, 

and 

[Ap, K), {K, Ap+i] if this last interval not empty. 

Proof. Let us find if there is an eigenvalue of J, that is a root of A^(A), in 
the interval [Xp,K). Using ([22]) and we get N{K) > > 0. From 
Lemma m we have two possibilities at Ap . Either G{Xp,n,n) = or 00. If 
the first happens, then Ap coincides with an eigenvalue of J by Lemma [5] 
(and N{Xp) = 6 hj fl36|l ). If the second possibility holds, then from (1321) 
A^(Ap+) = —00 . Since N is continuous in (Ap, K) there is at least a zero 
of in this interval. Therefore there is at least one eigenvalue of J in the 
interval [Ap, K). 

Now consider the interval [Ap_i,Ap). The two options mentioned above 
for Ap give us either iV(Ap) = 6' or A^(Ap— ) = 00. The two options for Ap_i 
are A^(Ap_i) = 9 and in this case Ap_i is eigenvalue of J, or iV(Ap_i+) = —00. 
If the second possibility happens, since A^(A) > for A near Ap, again from 
the continuity of A^(A) in (Ap_i, Ap) follows that has at least a zero in 
this interval. Therefore, there is at least one eigenvalue o f J in [Ap-_i,Ap). 
Continuing in this way, we get one eigenvalue of J in each interval [Xj, A^+i) 
with j = 1, ...,p — 1. Therefore if > Ap we get at least one eigenvalue of J 
in each [Xj, Aj+i) with j = 1, ...,p — 1 and one in [Ap, K), a total of at least 
p eigenvalues of J in the interval [Ai, K). 

Now consider the case K < Ap+i and let us see whether there is an 
eigenvalue of J in (K, Ap+i] . From ([32]) and ([37D , N{K) > 9. Now, from 
Lemma m we have two possibilities at Ap+i. Either G{Xp^i,n,n) = or 00. 
If the first happens, then Ap+i coincides with an eigenvalue of J by (!34l) 
(and A^(Ap+i) = 6* by (I5B]) ). If the second possibility holds, then from (15^ . 
A^(Ap+i— ) = —00 . Since is continuous in {K, Ap+i) there is at least a zero 
of N in this interval. Therefore there is at least one eigenvalue of J in the 
interval (i^, Ap+i]. Now consider the interval (Ap+i,Ap+2]. The two options 
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mentioned above for Ap+i give us either A^(Ap+i) = 6* or A^(Ap+i+) = oo. The 
two options for Ap+2 are A^(Ap+2) = ^ and in this case Ap+2 is eigenvalue of J, 
or A^(Ap+2— ) = —oo. If the second possibihty happens, since iV(A) > for A 
near Ap+i, again from the continuity of iV(A) in (Ap+i,Ap+2) follows that N 
has at least a zero in this interval. Therefore, there is at least one eigenvalue 
of J in (Ap+i,Ap+2]- Continuing in this way, we get one eigenvalue of J in 
each interval (Aj, A^+i] with j = p + 1, . . . , N — 1. 

Therefore if < Ap+i we get at least one eigenvalue of J in each (Aj, A^+i] 
with j = p + 1, . . . , N — 1 and one in {K, Xp+i], a total of at least N — p 
eigenvalues of J in the interval {K,Xn]- Since the p eigenvalues of J in 
[Ai, K) plus the N — p eigenvalues of J in {K, Xjy] give all the eigenvalues of 
J, we conclude that at most there is one eigenvalue of the perturbed operator 
in each one of the intervals considered. 

In case K = Ap+i we analyze first the interval {K, Ap+2] exactly as above, 
and find at least one eigenvalue of J in it. Continuing with the other intervals 
as before, we conclude that there is at least one eigenvalue of J in (Aj, A^+i] 
with j = p + 2, N — 1. Therefore we get N — p — 1 eigenvalues of J 
in {K,Xn]. These plus Ap+i and the p eigenvalues in [Ai,i^') give all the 
N eigenvalues of the perturbed operator. Therefore there is at most one 
eigenvalue of J in each of the intervals considered. □ 

Lemma 6. If K E <y{J) fl cr(J), then the following alternative holds: 
either 

(a) N{K) = e'^, and then N'{K) = 0, G{K,n,n) = and there are at 
most min(n — 1, — n — 2) other common points of cr(J) and o"(J). 

or 

(b) N{K) > 6"^, and then G{K,n,n) = 00 and there may be min(n, — 
n — 1) other common points of a{J) and cr(J). 

Proof. According to Lemma HI G{K,n,n) = or G{K,n,n) = 00. In the 
first case, by formula f l25|) . that means that the function 

1 9^ - N{X) 
X-K ' 1-^2 

has a zero at K, thus 

\{{K- A,) 

N{K) = = 

UiK-X,) 
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and, moreover, 



N'{K) = 0. 



Since G{X, n, n) may vanish only at min(n, N — n — 1) points of the spectrum 
of J, in the first case there may be at most min(n — l,N — n — 2) other points 
of a(J) na(J). 

In the second case, the function 9'^ — N{X) may not have a zero at K of 
order greater than 1, because otherwise G{X,n,n) would not have a pole at 
this point. We also know that G{X,n,n) has a negative residue at K. Thus, 
N{K) > 9^. In the second case there may be min(n, N — n — 1) more points 
of cr( J) n (t( J) where G(A, n, n) vanish. □ 

4 Inverse problem 

It turns out that the properties of the spectral data, described in the previous 
section, are sufficient. 
Let be given: 

\) a = {Ai, . . . ,Xn} and a = {Ai, . . . , Xn} two finite subsets of M Aj < 
Xi+i, Xi < Aj+i for i = 1, . . . , TV - 1. 
n) K eR 

iii) An integer number n E {0,1, . . . , N — 1} 

Introduce the following notations: 
{/xi, . . . , J := (7 n 0- - {K}. liana- {K} = 0, then g := ; (42) 



9^ ■= N{fii). liana-{K} = 0, then fix any p e {0,N{K)] and set 9^ := p. 



N{X) : 



nf=i(A-A,) 
n7=i(A-A,) 



(43) 



(44) 



n :— min(n, N — n — 1) 
Theorem 3. (Necessary and sufficient conditions) 



(45) 



The conditions 
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I) a, a and K interlace as in Theorem 

//; = iv(/ii) = iv(/i2) = ... = iv(^,) G (0, i) 

///; IfK^aUa then q < h and N{K) = 9^. 

IV) If K ^ <tU a then K E cr D a and either 

a) q < h and N{K) > 9"^ or 

b) q < n and N{K) = 9"^, N'{K) = (' denotes derivative) 

are necessary and sufficient for the existence of N x N Jacobi matrices J 
and J, where J is obtained by perturbing J at the n place as described in 
that IS an = 9'^{an + M), K.-i = 6'6„„i, 6„, = 9bn with M = {9^'^ - 1)K, 
such that 

a = spectrum of J, a = spectrum of J. 

Remark 2. Observe that condition I) implies N{K) G (0, 1) 

Theorem 4. Assume conditions of previous theorem hold. If K ^ a U a or 
K E a U a and option IV a) happens, then there are infinitely many pairs 
J, J of N X N Jacobi Matrices, if q ^ 0. Indeed this inverse spectral family 
is a collection of 

'N -2q- V 



n — q 

disjoint manifolds of dimension q and diffeomorphic to a q dimensional open 
ball. If q = 0, then there is only the finite 

iV- 1 
n 

number of pairs J, J. 

If K E a U a and option IV)b happens, there are infinitely many pairs. The 
inverse spectral family is a collection of 

N -2q-3 
n — q — 1 

disjoint manifolds of dimension q + 1 diffeomorphic to a q + 1 dimensional 
open ball. 
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Remark 3. Note that when K ^ a U a and in case IVb) the parameter 9"^ 
is uniquely determined by the spectral parameters i-iii), and in case IVa) 
9"^ is arbitrary in (0,A''(i^)). In particular, in Theorem H] this means that 
if cr U o" = {K} {q = 0) and N'{K) = 0, then there are a finite number 
(^-^) of solut ions for each 9'^ G {0,N{K)) and a collection of disjoint 

one-dimensional manifolds of solutions for 9^ = N{K). 

Proof. Here we prove simultaneously Theorems [3] and HI 

The necessity of the conditions I)-IV) is already proved in previous sec- 
tion: Theorem [2] proves necessity of condition I). Assertions (IMl) and fl36|) of 
Lemma O prove necessity of condition II). For condition III) use assertions 
fl55]) and fl5^ of Lemma O First part of condition IV) is fl3^ . Lemmata) 
implies condition IV)b. Lemma [6] b) and (l37j) . imply IV) a. Now we prove 
the sufficiency part of theorem [3] and theorem H] by finding all pairs of Jacobi 
matrices that have the given spectral data. 

Consider the function 

(compare to formula fl25|) ). Let us now prove that this is the Green's function 
of a Jacobi matrix. We consider two cases: 
Case A) K ^aUa 

Expanding in partial fractions we get: 
where 



/j,= lim^(A-A,)= nf=i(A.-A.) 

^ - n5,(A,-A.)(A,-ir) 



A^A, X - K 

and 



(3o = N{K) (49) 
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From condition III) and fH^ we obtain 

Now, from fH6|) and fl50l) we find that 
lim AG'(A) : 

A— j-oo 

respectively. Therefore 



and 

G(A) 

where := 

From f HHj) we know that /3j = if and only if Aj G crrio". From condition III) 
it follows that the sum in (!53ll has k := N — q terms and N — h < k < N . 
Using the expression (HHl) and the interlacing condition I), it follows that if 
ftj 7^ then > 0. Therefore 

fc fc 
G(A) = V with V = 1 and a^, > (54) 

1=1 ~ ^ 1=1 

According to theorem 6.2 of [3j, ( 15^ implies that G{X) has the form of a nn 
Green's function for at least one Jacobi matrix J. Now we shall describe the 
family of matrices which have G{X) as its nn Green's function and moreover 
have spectrum equal to the given set a. 

All finite Jacobi operators with a nn Green's function given by G{X) in 
f l5^ have the same eigenvalues X^,l = 1, ...,k, but the other N—k eigenvalues 
may change. To study the family of operators which correspond to a given 
Green'function, we shall use the theory of interior inverse problems for finite 
Jacobi matrices, developed in |3], theorems 6.1-4. The key formula of this 
method is (2.18) in [3]: 

- (J - zy^5ry^ = A - a„ + 6^m+(A,n) + bl_-^m^{X,n), (55) 



-1 and lim AG(A) = --^VA 

A^oo i — U 

1=1 



(51) 



1 ^ 

i=l 



(52) 



N N 

EaT^ with J]a. = l 



1=1 



(53) 



1=1 
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where 

m+(A,n) := {J[n+i,N-i] - 

m_(A,n) := (J[o,„-i] - 2)"^5„-i) 

are the so called m-Weyl functions. The matrices J[n+i,N-i] and J[o,n-i] were 
defined in Lemma[2l It happens that m_|_(A, n) determines J[„+i Ar_i] uniquely 
(see Remark H]), and has the form 

N-l~n 

m+(A,n)= J2 TTTx' > ^ (^6) 

i=l 

where ^^7^~"'7j = 1 and the are the eigenvalues of J[ri+i,Af-i]- Any sum 
of this form is legal for m^{X,n). Similarly m_{X,n) determines uniquely 
J[o,n-i] and has the form 

n 

m_(A,n) = V^^, ^i>0 (57) 

where Ylli=i^i — ^ 9i eigenvalues of J[o,n-i]- Any such sum is 

allowed for m_(A,n). 

The reconstruction procedure is as follows: 
Given G{X) as in f l5^ then 



fe— 1 
1=1 ^' 



ui — X 



where vi < V2 < ... < Vk-i are the zeros of G'(A). The numbers z/;,a and 
/3; > are determined by ctj, and Xi^ in the expression Now we have to 
write the right side of equality (l58l) in the form of the right side of equality 
([55]) for some (3„,6^,6^_]^ and ??7,4.(A, Ti), 'n2_(A, Ti) of the form described in 
(!56|) and (!57|) . If we do this, then we would have according to (l55ll that G{X) 
is a nn Green's function for a matrix J with corresponding entries a„, hn-i 
and submatrices J[„+i_Ar_i] and J[o,n-i] determined by the m-Weyl functions. 
From condition II) we know that the q points of a fl a are among zeros of 
G{X). We will construct m+(A, n) and m_(A, n) in such a way these g points 
are common poles of them. The other k — 1 — q zeros of G{X) will be poles 
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of just one of the m-Weyl functions. Since m_(A,n) has n poles, then there 
are 

'k-l-q\ _ /N -2q-r 
n — q J \ n — q 

possibihties of distributing non common poles. For each of the q common 
ole 

o(i) 



(59) 



poles jjLi we pick a decomposition f3i = P^^^ + P^'^^uch that the addends 



= 1,2 appear each in one of the sums fl56|) . fl57j) . So we have q pa- 
rameters which generate a manifold for each one of the possible choices (159|) . 
That these manifolds are diffeomorphic to a sphere follows from theorem 
3.6 of [3]. Since we have constructed the m-Weyl functions, the matrices 
J[,„_i_i jv_i] and J[o,n-i] are determined. We can fix a = a„ and 



I so fj.j^ is anfj^ I so /j,^ is anfj^ I so /j,^ is an I so /j,^ is anfj^ 

fij^ is not an and an /j,^ is not an and an g^ 

where /j and Qi are defined in ( ISBj) . (IFT]) . So given G'(A) we have constructed 
a family of Jacobi matrices so that each member of it, has G{X) as its nn 
Green's function and its eigenvalues are exactly the points of a. We have only 
to prove that if we perturb one of these matrices J, then the perturbation J 
has spectrum exactly a. Consider then the Jacobi matrix J, obtained from 
J by formulas with M := (6'~^ — 1)K, where K and 6^ are as defined 
above in theorems Ei Then, by formulas fl2^ and 



1 r 9^ N{X) 1 1 f 9^ N{X) 



l-9nX-K X-K) l-9nX-K X-K. 
which implies N{X) = N{X) and n(-^ - ^i) = IK-^ - h)^ thus 

cr = {Ai, . . . , Aat} = {Ai, . . . , Xn} = cr{J). 
This proves all assertions of theorems [3] and H] for case A)K ^ a U a. 

Case B)K eaUa. 
In this case from condition IV) we know that there exist jo ^ {l? ^} such 
that K = Xjg = AjQ. Then (H7|l takes the form: 
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where 
and 

= NiK) (62) 
Analogously to what was done in fl^ . flS2]) and fl5^ we get 

1 ^ ^2 

T3^Eft-T3^-i (63) 

i=l 

and 

AT TV 

^(^) = Ea^ 5^«. = 1 (64) 

2=1 ^ i=l 

where at := ii i ^ jo ■, o.jo = ^'^^2 • For i ^ jo we get using fl6T|) that 
Qij = if and only if Aj = Aj, that is exactly q times according to fH2]) . Using 
the interlacing condition I) and fl6T!) it follows that a, > if 7^ 0. Therefore 

fc-l k-l 

G'(A) = 5^ + with ^ a, + = a,>0 (65) 
«=i 1=1 

Where k := N — q. Now we have two options: 

i) If situation IV) a happens, then we get exactly (15^ since ajg > 0, is a 
pole of G{X) and the analysis is completely analogous to Case A). 
i))If IV)b holds then aj^ = , then 

G{X) = V with ai, > (66) 

1=1 " ^ 

and iiT is a zero of G{X). So we have g + 1 fixed zeros. An analysis similar 
to the one for Case A) gives 

k-2-{q+l)\ _ fN -2q-3 
n — {q + 1) J \n — q — l 

possible choices and then g + 1 parameters. 

□ 
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Remark 4. To find the entries of the matrices J[o,n-i] and J[„+i,Ar-i] one 
could use the continuous fraction expansions of the Weyl functions: 

- m_(A,n)"^ = A - a„_i ^ _^ (67) 

hi 



\ n— 3 

A — a„_2 — 



A - Ctn-S- ■ • k 



A - ao 



and 



52 

m+(A, n)"^ = A - a„+i "'^^ ,0 (68) 

0: 



A — a„+2 



A - fln+s- ■ • b 



,2 

Ar-2 



A — CLj^—l 



These expansions are unique and can be obtained using Euchd's algo- 
rithm. We have to choose the negative square root of the since the off 
diagonal terms of our Jacobi matrices are negative, see formula ([2]). 
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